Topological phase in one-dimensional interacting fermion system 
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We study a one-dimensional interacting topological model by means of exact diagonalization 
method. The topological properties are firstly examined with the existence of the edge states at 
half-filling. We find that the topological phases are not only robust to small repulsive interactions 
but also are stabilized by small attractive interactions, and also finite repulsive interaction can drive 
a topological non-trivial phase into a trivial one while the attractive interaction can drive a trivial 
phase into a non-trivial one. Next we calculate the Berry phase and parity of the bulk system 
and find that they are equivalent in characterizing the topological phases. With them we obtain 
the critical interaction strengths and construct part of the phase diagram in the parameters space. 
Finally we discuss the effective Hamiltonian at large-[7 limit and provide additional understanding of 
the numerical results. Our these results could be realized experimentally using cold atoms trapped 
in the ID optical lattice. 



I. INTRODUCTION 

The finding of time-reversal invariant topological insu- 
lators (TIs) has become an exciting event in condensed 
matter physics. Since then lots of works have been car- 
ried out theoretically and experimentally, predicting and 
verifying many exotic physical properties exhibited by 
TIsCuEIl. The key feature of TIs is the existence of the 
robust edge states determined by the bulk topological 
property, which can be described by Zi— valued topo- 
logical invariants"^. The original definition of TIs is for 
non-interacting band structures and the relevant physics 
has been well understood. Thereafter one of the subjects 
that need to be explored further is the effects of electron 
correlations on TIs. 

In the situation where many-body interactions exist, 
the definition of TIs from the topological field theory, 
which is the presence or absence of a topological term in 
the effective electromagnetic action, is generally valid 8 . 
In addition, the method using the Green's functions to 
construct the topological invariants is also appliable^EE] 
Yet the main difficulty is still to deal with the interac- 
tions properly. At the mean-field level, it has been shown 
that interactions can change the trivial insulators into 
non-trivial ones^"^. It is also been proposed that TIs 
and new topological phases may appear in the systems 
with consid erable interactions such as 4c? or 5d transition 
metal oxides^UJil Recently several numerical simulations 
and analytical works were performed on the interacting 
two-dimensional Kane-Mele model and studied the inter- 
play of spin-orbit coupling and Coulomb repulsiorP^HUl. 
These studies show that the Hubbard repulsive interac- 
tion can transform the TI of Kane-Mele model to ei- 
ther the spin liquid phase or antiferromagnetic insulating 
phase depending on the strength of the spin-orbit cou- 
pling. A study using the Lanczos algorithm concludes 
that the topological properties have already manifested 
themselves in small systems and therefore can be studied 
numerically via exact diagonalization (ED) and observed 
exp eriment all jEJ 



There also appear works addressing the question that 
how the presence of interactions changes the classification 
of the topological phases. For noninteracting systems, 
five symmetry classes are topological nontrivial in each 
spatial dimensionality. This classification has been ex- 
pected to be also appliable to interacting systems as long 
as the strength of the i ntera ctions is sufficiently small 
as compared to the garj 26 * 27 [ However recent studies 
on a specific one-dimensional (ID) model show that the 
free-fermion classification breaks down in the presence of 
interaction s] 28 * 29 ^. So it is most possible that the inter- 
action doesn't modify the topological nontrivial classes 
uniformly. 

In this paper, we use the ED to study the effect of in- 
teractions in a ID lattice model which is known to have a 
topological state in its free form. With the ground-state 
energies and wave functions, we firstly calculate the en- 
ergy and the distribution of the quasi-particle added or 
removed from the system at half-filling for open (OBC) 
and periodic boundary conditions (PBC), then identify 
the topological phases with the existence of the edge 
states. The topological features indeed have manifested 
themselves clearly in the small sizes we can access. We 
consider both the repulsive and attractive interactions 
and find the topological phase is robust to repulsive in- 
teraction while stabilized by attractive one when the in- 
teraction strengths are small. For repulsive interaction 
we find that as its strength is increased, the system un- 
dergoes a topological quantum phase transition (TQPT) 
into a trivial insulator. For attractive interaction it can 
drive a trivial insulator into a non-trivial one. Then 
we calculate the Berry phase and the parity of the bulk 
system, which are equivalent in characterizing the topo- 
logical property of the system. When the system is in 
the topological phase, they have nontrivial values. With 
them we obtain the critical interaction strengths and con- 
struct part of the phase diagram in the parameters space. 
Finally we discuss the effective Hamiltonian at large-C/ 
limit and provide additional understanding of the numer- 
ical results. 
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FIG. 1: (Color online) (a)The tight-binding band structure 
of H(k). The low-energy excitation is located at k\ (black 
line) or fej (green line) depending on the parameters, (b) 
Edge modes in the topological phase (M = -1) on a chain 
of length N — 20 with OBC. (c) Berry phase of the occupied 
Bloch state of Eq.{T| at different M. (d) The plot of the curve 
[d x {k),d z (k)] with k £ [0, 2tt]. In all figures B = 1 is set. 



II. THE ID MODEL WITH EDGE MODES 

Our starting point is the ID non-interacting tight- 
binding modelj^QI 

H = £(M + 2fl)*Jff,tf i -£B*J<7 J( ¥ i+4 (1) 

i i,x 

- sgn(x)iA&l(Tx i &i+i 

where <r x , a z are Pauli matrices and ty{ = (ctf, c^) T with 
c it( c «i) electron annihilating operator at the site r^. The 
first two terms represent the differences of the on-site 
potentials and the hopping amplitudes between the up- 
and down- electrons, and the third term is due to the 
spin-orbit coupling. In momentum space Eq.([T]) becomes 

#o = with *fe = {ckt,Cki) T the Fourier 

partner of ^ and 

H(k) = [M + 2B - 2Bcos(k)]<r z + 2Asin(k)a x . 

The spectrum of H(k) consists of two bands, 

E£' 2) = ±y/[M + 2B- 2Bcos(k)] 2 + [2Asin(k)} 2 . 

Usually the two bands are dispersive, but when the pa- 
rameters satisfy —M/2 = B = A the bands are flat. For 
M = (M = -AB) bands 1 and 2 touch at the Dirac 
point k\ — (&2 = tt), while for other values a gap 



A=min{2|M|, 2\M + AB\} opens up at the Dirac point 
k% or k2- At half- filling, depending on the values of the 
parameters A, B and M the system can be a trivial insu- 
lator or a non-trivial insulator with edge modes. In the 
following of the paper we take B positive and set A = 1 
as the energy scale. 

The topological property of the system can be under- 
stood in terms of Berry phase in k space, which is 7 = 
§ A(k)dk with the Berry connection A(k) — i{uk\-^\uk) 
and \v,k) the occupied Bloch stat d 32 * 33 -! The numerical 
result is shown in Figfljc). It shows the Berry phase 7 
mod 2-7T gets a nonzero value ir for — AB < M < 0. We 
also have performed numerical diagonalization of H$ on 
a chain with OBC. In accord with the above argument, 
we find a pair of zero modes appearing in the gap when 
the Berry phase of the system is 7r (Figfljb)). 

We notice that when the Berry phase is tt the masses 
at the two Dirac points fci and /C2 have different signs, 
i.e., M(M + AB) < 0. This can also serve as a cri- 
terion of the topological property in the system. The 
reason can be understood from the low-energy Hamil- 
tonians gov erning the excitations in the vicinity of the 
Dirac points^ 31 * 36 * 37 *. By linearizing H(k) near k\ and 
ki we obtain two Dirac Hamiltonians, 



hi = 2Aa x k + 
hi = -2Aa x k 



+ (Ar- 



ea) 



AB)<j 2 



For — AB < M < the masses at the Dirac points ex- 
hibit opposite signs, while for M > and M < —AB the 
masses at the Dirac points exhibit the same signs. In the 
following we show that the two cases correspond to two 
phases with different topological properties. We consider 
a junction between the two phases running along a line in 
real space (suppose — AB < M < for x < and AT > 
for x > 0). So the mass of h\ necessarily undergoes a 
sign change across the x = boundary. Such a soliton 
mass profile is known to produce massless state in the 
associated Dirac equation, localized near the boundary. 
Specifically, Dirac equation 

[2A{-i)a x d x + cr z m,(x)]<p(x) = E<j>{x) 

with m(x — > —00) < and m[x — > 00) > has a gapless 
solution 



\ — f m(x')dx' , 1 \ 

z) — e 2A Jo K ' ( . ) 



localized at the boundary with zero energy. Since the 
M > phase can be continually connected to the M = 00 
phase which is a trivial insulator, the phase for —AB < 
AT < has non-trivial topological property. The above 
argument is similar when M < —AB for x > in the junc- 
tion. So through the relative signs at the Dirac points we 
can obtain the same condition for the non-trivial topo- 
logical phase. 

Furthermore, we can write H(k) = d(fc) • <?, where 
d x (k) = 2Asin(k) ,d y (k) = and d z (k) = M + 2B - 
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FIG. 2: (Color online) The energies of the quasi-particles added or removed from the half-filling system with OBC and PBC. 
(a) and (b): the phase at U = is topological (M = —1); (c): the phase at U — is trivial (M — 1). The critical interactions 
U c in (a) and (c) can be accurately determined from the parity of the ground-state wave function. Here B = 1. 



2Bcos(k). The pair [d x (k),d z (k)] forms a closed loop in 
the plane when k changes from to 2ir. If the system 
has non-trivial topological properties, the loop will con- 
tain the origin point of the plane (Figjljd)). It is also 
consistent with the former arguments. 



III. THE EFFECT OF HUBBARD 
INTERACTION 

To study the effect of interaction, we add the Hubbard 
term Hi = UJ2i n it n 4 *° Eq.Q. This term changes the 
single electron Hamiltonian EqTfl]) to a many-body one. 
For non-interacting case, the topological property can 
be directly read from the presence of the edge states in 
the gap of the single electron energy spectrum. However 
in the presence of interaction, there is no longer single 
electron state. Then how can we identify the topological 
property of the system? 



A. The existence of the edge states 

Generally for a finite chain of N sites, the full Hilbert 
space of the system has dimension A N . Since in the 
Hamiltonian Hq + Hi the total particle number is con- 
served, we can get the ground-state energy and wave 
function of a system with a fixed number of electrons 
using ED. Then we can define the energy of the quasi- 
particle added to a system with n electrons as AE n = 
E ( r \ +1 — E®, where E^ is the ground energy of a system 
with n particles. Similarly to the non-interacting case, if 
the system has non-trivial topological property, there ap- 
pear states in the gap of the quasi-particle energy spec- 
trum (QPES) as the boundary condition changes from 
PBC to OBC. Since this definition can be continually 
connected to the single electron case, we expect it valid 
at least for small-l/ cases. 



Firstly we consider the effect of repulsive interaction. 
Figure [2]ja) shows the results of starting from a system 
with non-trivial topological property. Similar to the non- 
interacting case, we concern about the electron added or 
removed from the half-filling system. It shows that the 
energies of the quasi-particles added or removed appear 
in the gap when U is below a critical value U c . At small U 
the two in-gap modes have exactly the same values. Due 
to the finite-size effect, when U approaches the critical 
value, the energies of the two modes is separated by a 
gap. When U is beyond U c , the in-gap modes disappear 
and evolve into the bulk ones. The result clearly shows 
that the topological phase survives in the presence of 
small repulsive interaction and a TQPT is driven by finite 
repulsive interaction. We also study the case of starting 
from a trivial insulator and find that no in-gap modes 
appear when repulsive interaction is added. 

Next we turn to study the effect of attractive inter- 
action. The calculations are straightforward and the re- 
sults are shown in Figj2]jb) and (c) . In Fig[2^b) attractive 
interaction is introduced into a system with non-trivial 
topological property. It shows that the existence of the 
in-gap modes persist to quite large strengths. Then as 
the strength is further increased, the in-gap modes con- 
tinuously evolve into bulk ones. In Figj2jc) we show the 
result of introducing attractive interaction to a trivial in- 
sulator. We find that when the strength reaches a critical 
value U c , the in-gap modes begin to appear in the gap 
of QPES, indicating that attractive interaction can drive 
a trivial phase into a non-trivial one. Then the in-gap 
modes persist till the strength becomes very large when 
they evolve into b ulk o nes. This behavior is very similar 
to that of disordeiPS. 

Till now by identifying the in-gap states in the QPES, 
we show the effects of repulsive and attractive Hub- 
bard interactions in the topological phase. In the non- 
interacting systems, the in-gap mode is also referred to 
edge mode due to the fact they mainly distribute near 
the edges of the chain. Similarly in the interacting case 
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FIG. 3: (Color online) The distribution of the quasi-particles 
added or removed from the half-filling system with OBC. (a) 
and (b): M = —I; (c) M = 1. Here B = I and the system 
size N = 10. 



we can also calculate the distribution of the in-gap mode 
to study its nature. Using the many-body wave func- 
tions, the distribution of the electron added can be de- 
fined as: An, = ~ (V>£l™ilV>£}> where 
hi = c].j.Cj<f + ct^Cjj, is the electron number operator on 
site i and i/ 1 ^ is the ground-state wave function of the 
system with n electrons. The results at different U and 
M are shown in Fig[3j Since the electrons added or re- 
moved at half-filling have exactly the same distributions, 
only one of them is shown. In Figj3ja), we start from 
a system with nontrivial topological property. As ex- 
pected at U — the in-gap state mainly distributes near 
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FIG. 4: (Color online) The energies and parities of the 
ground- and the first-excited states. (a)M = —1 and (b)M = 
1. The red (black) curve has the parity value -1(1). In both 
figures B = 1 and N — 8. 



the edges. As we increase the strength of the repulsive 
interaction, the distribution begins to evolve from the 
ends to the bulk. Though the result is greatly affected 
by the finite-size effect, the distribution at U < U c is 
still clearly distinct from that at U > U c (U c is about 
3 ~ 4), indicating the different topological properties ex- 
isting in the system. Next we add attractive interaction 
to the above system and the result is shown in Fig[3^b). 
As the strength increases, the topological phase is firstly 
stabilized, manifested by the increase of the components 
near the ends. Then the distribution begins to oscillate 
between positive and negative values. The occurrence 
of negative value is due to the many-body effect. The 
amplitude of the oscillation decays from the end to the 
bulk. And as the strength is further increased, the de- 
cay becomes slower, so the in-gap state begins to disap- 
pear, indicating the system experiencing a TQPT. We 
also start from a trivial insulator and the result is shown 
in Figj3]jc). It shows that in the presence of repulsive 
interaction the distribution still mainly concentrates in 
the bulk, while after attractive interaction is added the 
distribution begins to evolve from the bulk to the ends. 
The dramatic change occurs at U = —3 when the distri- 
bution shows a sign of edge one. The estimated critical 
value —3 ~ —2 is in good consistent with U c in Figj2|c). 
Then as the strength of the attractive interaction is fur- 
ther increased the distribution shows a similar behavior 
as that in Figj3^b) . Our these results are in good consis- 
tent with those obtained from the QPES. 



B. The topological properties of the bulk system 

Till now we have examined the topological phase in the 
interacting system with the existence of the edge states. 
It is known that their existence is due to the bulk topolog- 
ical properties. So in the following we will calculate the 
Berry phase of the ground-state of the interacting system 
at half-filling using the twisted boundary conditionsP^IUl. 
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FIG. 5: (Color online) The phase diagram in (M, U) and 
(M, B) planes. The system size N — 8. 

It can be denned as 

where 9 is the twisted boundary phase which takes val- 
ues from to 2ir and ipg is the corresponding ground- 
state many-body wave function at half-filling. To com- 
pare with the previous results, we firstly study the Berry 
phases associated with Fig(2j We find that the Berry 
phase 7 mod 2ir has the value it and for U < U c and 
U > U c (see Figg}. 

The Hamiltonian %{k) possesses the symmetry imple- 
mented by a unitary transformation 



a z H{k)cr z = H(—k). 



(3) 



The symmetry is similar to the inversion symmetry ex- 
cept that there is an additional sign when it is performed 
on the spin-down electron. Thus we can define the par- 
ity of the wave function ijQ. In the basis of fixed num- 
ber of electrons, tp® — X^<^N) wrtn I*) denoting the i- 
th n-electrons basis. Under the inversion transforma- 
tion, Tip° n = Y l i<t>i'P\ i ) = Y,i<t>iZi\j)i = ±^°) where 

V\i) = and e$ = ±1 depending on the times of 

exchanging the fermion operators. We calculate the par- 
ities of the ground- and the first-excited states associated 



with Fig(2] and the results are shown in Figj4] It shows 
that at half-filling there is a gap between the two states 
even in the presence of interactions and the gap may be 
eliminated at the critical strength U c where the parity 
also change its sign. 

We also find that the parity is the same as the Berry 
phase to characterize the bulk topological properties 
since they both change their values at the same critical 
value U c . From either one the accurate values of the crit- 
ical interactions can be obtained, such as: U c = 3.66 in 
Fig{2ja) and U c = —2.86 in Fig{2jc). We want to mention 
that in some cases such as at large attractive interaction 
though the in-gap states disappear, the Berry phase or 
the parity of the ground-state doesn't change their val- 
ues. These cases are beyond our discussion because the 
ground-states become degenerate and the systems aren't 
insulators any more. 

From these calculations, we can get part of the phase 
diagram in the parameters space, as shown in Fig(5] 
Fig(5ja) is in (M, 17) plane, it shows that at each value 
of B the curve is symmetric about M = —2B and the 
behaviors are similar as M is far from — 2B. When M 
approaches —2B, the property is different depending on 
the value B. For B < 1, the critical interaction gets its 
maximum at M = —2B. While for B > 1 the results 
become complex and depend greatly on the system sizes, 
so no definite conclusions are made. Fig[5jb) is in (M, B) 
plane. At U = the topological phase is in the area re- 
stricted by the lines M = and B = — M/4. Its size is 
shrunk by repulsive interaction and broadened by attrac- 
tive one. Corresponding to the results in Figj5](a) at fixed 
B of each curve the sum of the two corresponding M is 
— 4B and when U > and B < 1 the tips of the curves fall 
on the line B = -M/2. We notice that (M, B) = (-2, 1) 
is a special point where the non-interacting system ex- 
hibits flat bands. Here at large repulsive interaction the 
ground- and the first-excited states also become degen- 
erate and the Berry phase and parity of the ground-state 
keep the same value all the way. 



C. The effective Hamiltonian at large-[/ limit 

To understand the phase diagram better, it is helpful 
to study the system at large-[7 limit, when the effective 
Hamiltonian ia^l, 



H, 



eff 



E 



(Jxcr, <J l+1 



J 2 « +1 )(4) 



Y J {M + 2B)al 



where J x = (B 2 - A 2 )/U,J V = (B 2 + A 2 )/U and 
Jz = - Jx- When A = B it is H l ef} = ./ V,a;V f . • + 

Jg of with J = 'jj- and g — ^ M ^gf^ U , which is 
the quantum Ising modeP^. By tuning the dimension- 
less coupling g, i?e// exhibits a quantum phase transi- 
tion and the critical point is exactly at \g\ = 1. For 



G 



2B < M < ±fj 2B, the ground state is 

twofold degenerate and possess long-range correlations in 
the magnetic order parameter a y , while beyond the above 
range of M it is in a quantum paramagnetic state. This 
naturally explains the degeneracy at B = — M/2 = 1 for 
large repulsive interaction. 

IV. CONCLUSIONS 

To conclude, we study the effect of interactions in a 
ID topological model by means of the ED method. The 
topological features have already manifested themselves 
clearly in the small sizes the calculations can access. Our 
studies focus on the half-filling system. We examine the 
topological phases with the existence of the edge states, 
which arc exhibited from the energy and distribution of 
the electron added or removed to a system at half-filling 
and with OBC. We show that the topological phase is ro- 
bust to small interactions and finite repulsive interaction 
can drive a topological non-trivial phase into a trivial one 
while the attractive interaction can drive a trivial phase 
into a non-trivial one. 

We calculate the Berry phase and the parity of the 
ground-state wave function to study the bulk topological 
properties. The Berry phase and the parity have intrinsic 
connection and are equivalent to describe the topologi- 
cal properties due to the symmetry in our model. In the 
cases where the edge state exists, the Berry phase has 
nontrivial value ir. At the TQPT points the Berry phase 



and the parity change their values. From them we deter- 
mine the critical interactions and construct part of the 
phase diagrams in the parameters space. 

Our these results demonstrate the existence of the 
topological phases in ID interacting fermion systems. 
Though the model we use is artificial, it may be con- 
structed experimentally using cold atoms trapped in the 
optical lattice, which allows one to directly simulate ideal 
and tunable models. At present ultracold Fermi gases in 
a truly ID regime can be realized using strong optical lat- 
tices and the interactions can be tuned between repulsive 
and attractive ones by means of Feshbach resonances^. 
With these developments some of the basic phenom- 
ena in the interacting fermion systems are being stud- 
ied. Besides, there appear some studies which suggest 
the methods to mimics the effect of spin-orbit coupling 
and produce topological states of matter in cold-atom 
systems^H^l So it is very hopeful that our these results 
are tested in cold- atom experiments. 
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